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Abstract. For continuous-variable systems, we introduce a measure of entanglement, the 
continuous variable tangle (contangle), with the purpose of quantifying the distributed (shared) 
entanglement in multimode, multipartite Gaussian states. This is achieved by a proper convex 
roof extension of the squared logarithmic negativity. We prove that the contangle satisfies 
the Coffman-Kundu-Wootters monogamy inequality in all three-mode Gaussian states, and 
in all fully symmetric TV-mode Gaussian states, for arbitrary N . For three-mode pure states 
we prove that the residual entanglement is a genuine tripartite entanglement monotone under 
Gaussian local operations and classical communication. We show that pure, symmetric three- 
mode Gaussian states allow a promiscuous entanglement sharing, having both maximum 
tripartite residual entanglement and maximum couplewise entanglement between any pair of 
modes. These states are thus simultaneous continuous-variable analogs of both the GHZ and 
the W states of three qubits: in continuous-variable systems monogamy does not prevent 
promiscuity, and the inequivalence between different classes of maximally entangled states, 
holding for systems of three or more qubits, is removed. 
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One of the main challenges in fundamental quantum theory as well as in quantum 
information and computation sciences lies in the characterization and quantification of 
bipartite entanglement for mixed states, and in the definition and interpretation of multipartite 
entanglement both for pure states and in the presence of mixedness [l] El- More intriguingly, 
a quantitative, physically significant, characterization of the entanglement of states shared 
by many parties can be attempted: this approach, introduced in a seminal paper by 
Coffman, Kundu and Wootters (CKW) 1 3 1, has lead to the discovery of so-called "monogamy 
inequalities", constraining the maximal entanglement distributed among different internal 
partitions of a multiparty system. Such inequalities are uprising as one of the fundamental 
guidelines on which proper multipartite entanglement measures have to be built |4|. 

While important insights have been gained on these issues in the context of qubit systems, 
a less satisfactory understanding has been achieved until recent times on higher-dimensional 
systems, as the structure of entangled states in Hilbert spaces of high dimensionality exhibits 
a formidable degree of complexity. However, and quite remarkably, in infinite-dimensional 
Hilbert spaces of continuous-variable systems, important progresses have been obtained in the 
understanding of the entanglement properties of a restricted but fundamental class of states, 
the so-called Gaussian states 1 5 6 1 . These states, besides being of great importance both from 
a theoretical point of view and in practical applications, share peculiar features that make their 
structural properties amenable to accurate and detailed theoretical analysis Q. 

In this work we address the problem of distributing entanglement among multiple 
modes of a continuous variable system. We introduce the continuous-variable tangle to 
quantify entanglement sharing in Gaussian states and we prove that it satisfies the Coffman- 
Kundu- Wootters monogamy inequality 0. Nevertheless, even in the basic instance of three 
modes, we show that pure, symmetric Gaussian states, at variance with their discrete-variable 
counterparts, allow a promiscuous sharing of quantum correlations, exhibiting both maximum 
tripartite residual entanglement and maximum couplewise entanglement between any pair of 
modes. 

The paper is organized as follows: in Sec. ^ we review the basic properties of 
Gaussian states of continuous variable systems, and set up notations; in Sec.|2]we address 
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the quantification of entanglement sharing in such states, introducing a new entanglement 
monotone which is shown to generalize the tangle defined in discrete-variable systems; 
in Sec. [3] we apply this measure to prove that all three-mode Gaussian states and all 
symmetric multimode Gaussian states satisfy a monogamy inequality for continuous-variable 
entanglement, and that in the specific case of three-mode states the residual entanglement, 
emerging from the monogamy inequality, is a genuine tripartite entanglement monotone; in 
Sec. |4] we exploit this result to investigate the sharing structure of tripartite entanglement 
in Gaussian states, unveiling striking differences with their discrete-variable counterparts; 
finally, in Sec. [5] we summarize our results and outline possible roadmaps ahead. 

1. Gaussian states: structural properties 

In a continuous variable (CV) system consisting of N canonical bosonic modes, associated to 
an infinite-dimensional Hilbert space, and described by the vector X of the field quadrature 
operators, Gaussian states (such as coherent, squeezed, thermal, and squeezed thermal 
states) are those states characterized by first and second statistical moments of the canonical 
quadrature operators. When addressing physical properties, like entanglement, that must be 
invariant under local unitary operations, first moments can be neglected and Gaussian states 
can then be fully described by the 2N x 2N real covariance matrix (CM) cr, whose entries 
are o~ij = l/2({jQ, Xj}} — (Xi)(Xj). This allows, for Gaussian states, to indicate them 
indefferently by the density matrix p or by the CM er. A physical CM er must fulfill the 
uncertainty relation 

(r + m>o, (i) 

with the symplectic form fl = ©™ =1 w and ui = <%-i — <%+i, i,j = 1,2. Symplectic 
operations (i.e. belonging to the group Sp^N,]R) = {$ <= SL(2N, R) : S T QS = ft}) acting 
by congruence on CMs in phase space, amount to unitary operations on density matrices in 
Hilbert space. In phase space, any TV-mode Gaussian state can be written as er = S T uS, 
with v = diag {m, ni, n2, n2, . . . , rijv, tin}. The set £ = {rii} constitutes the symplectic 
spectrum of cr and its elements must fulfill the conditions rii > 1, ensuring positivity of the 
density matrix p associated to er. The symplectic eigenvalues rii can be computed as the 
eigenvalues of the matrix |if2<x|. The degree of purity p = Tr p 2 of a Gaussian state with CM 
er is simply p = 1/ \J Dct cr. 

Concerning the entanglement, positivity of the partially transposed state p (from now on 
"~" will denote partial transposition), obtained by transposing the reduced state of only one 
of the subsystems, is a necessary and sufficient condition (PPT criterion) of separability for 
(N + l)-mode Gaussian states of (1 x 7V)-mode bipartitions H]|9) and for (M + iV)-mode 
bisymmetric Gaussian states of (M x iV)-mode bipartitions 1101 . In phase space, partial 
transposition in a (1 x 7V)-mode bipartition amounts to a mirror reflection of one quadrature 
associated to the single-mode party If {hi} is the symplectic spectrum of the partially 
transposed CM cr, then a (N + l)-mode Gaussian state with CM cr is separable if and only 
if hi > 1 Vi. This implies that a proper measure of CV entanglement is the logarithmic 
negativity II II 

£?AT = In||p||i, (2) 

where j| • j| i denotes the trace norm 1121 . The logarithmic negativity £V is readily computed 
in terms of the symplectic spectrum hi of er as 

E M = - J2 Inn*- (3) 

i:hi <1 
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Such a measure quantifies the extent to which the PPT condition is violated. For two-mode 
symmetric states, the logarithmic negativity is equivalent to the entanglement of formation 
(EoF) E F in): 

E F (p) = inf Vp,B(W), (4) 
{pi,W'i}} '-r 

where E(\ipi)) is the von Neumann entropy (or entropy of entanglement) of the pure state \ipi), 
and the infimum is taken over all possible pure states decompositions p = ^ pi\ipi) (ipi\- 

In fact, the logarithmic negativity is positive defined, additive, monotone under local 
operations and classical communication (LOCC) II14I . constitutes an upper bound to the 
distillable entanglement in a quantum state p, and is related to the entanglement cost under 
PPT preserving operations 1151 . 

2. The contangle 

Our aim is to analyze the distribution of entanglement between different (partitions of) modes 
in CV systems. In Ref. Coffman, Kundu and Wootters (CKW) proved for system of three 
qubits, and conjectured for N qubits (this conjecture has now been proven by Osborne and 
Verstraete 1161 ), that the bipartite entanglement E (properly quantified) between, say, qubit 
A and the remaining two-qubits partition (BC) is never smaller than the sum of the A|B and 
A|C bipartite entanglements in the reduced states: 

E A\{BC) > E A\B +E A\C_ (5) 

This statement quantifies the so-called monogamy of quantum entanglement 1171 . in 
opposition to the classical correlations, which are not constrained and can be freely shared. 
One would expect a similar inequality to hold for three-mode Gaussian states, namely 

E i\(jk) _ E i\j _ E i\k > j (6) 

where E is a proper measure of bipartite CV entanglement and the indexes {i, j, k} label the 
three modes. However, the demonstration of such a property is plagued by subtle difficulties. 
Let us for instance consider the simplest conceivable instance of a pure three-mode Gaussian 
state completely invariant under mode permutations. These pure Gaussian states are named 
fully symmetric, and their standard form CM (see Refs. I18II19I ). for any number of modes, is 
only parametrized by the local mixedness ai oc = 1/ [ii oc , an increasing function of the single- 
mode squeezing r; oc , and ai oc — ► 1 + when r\ oc — > + . For these states the inequality 
can be violated for small values of the local squeezing factor, using either the logarithmic 
negativity E^r or the EoF Ep (which is computable in this case, because the two-mode 
reduced mixed states of a pure symmetric three-mode Gaussian states are again symmetric) 
to quantify the bipartite entanglement. This fact implies that none of these two measures 
is the proper candidate for approaching the task of quantifying entanglement sharing in CV 
systems. This situation is reminiscent of the case of qubit systems, for which the CKW 
inequality holds using the tangle r, defined as the square of the concurrence j20;|, but can fail 
if one chooses equivalent measures of bipartite entanglement such as the concurrence itself or 
the entanglement of formation Q . 

It is then necessary to define a proper measure of CV entanglement that specifically 
quantifies entanglement sharing according to a monogamy inequality of the form (jS). A 
first important hint toward this goal comes by observing that, when dealing with 1 x JV 
partitions of fully symmetric multimode pure Gaussian states together with their lxl 
reduced partitions, the desired measure should be a monotonically decreasing function / 
of the smallest symplectic eigenvalue n_ of the corresponding partially transposed CM er. 
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This requirement stems from the fact that n_ is the only eigenvalue that can be smaller 
than 1 1181 . violating the PPT criterion with respect to the selected bipartition. Moreover, 
for a pure symmetric three-mode Gaussian state, it is necessary to require that the bipartite 
entanglements E l I to'*) and-E' 1 '-' = E^ k be respectively functions f(n_ ) and/(n!i"') of the 

associated smallest symplectic eigenvalues fi 1 ^ 1 ^ and h l _[\ in such a way that they become 
infinitesimal of the same order in the limit of vanishing local squeezing, together with their 
first derivatives: 

f(h^ k) )/2f(^) ~ /'(n^ fe) )/2/'(^')-l for a loc - 1+ , (7) 

where the prime denotes differentiation with respect to the single-mode mixedness ai oc . The 
violation of the sharing inequality (|6j exhibited by the logarithmic negativity can be in fact 
traced back to the divergence of its first derivative in the limit of vanishing squeezing. The 
above condition formalizes the physical requirement that in a symmetric state the quantum 
correlations should appear smoothly and be distributed uniformly among all the three modes. 
One can then see that the unknown function / exhibiting the desired property is simply the 
squared logarithmic negativity 112 II 

/(n_) = [-lnn_] 2 . (8) 

We remind again that for fully symmetric (JV + l)-mode pure Gaussian states, the partially 
transposed CM with respect to any 1 x JV bipartition, or with respect to any reduced lxl 
bipartition, has only one symplectic eigenvalue that can drop below 1 1181 : hence the simple 
form of the logarithmic negativity (and, equivalently, of its square) in Eq. i|8}- 

Equipped with this finding, one can give a formal definition of a bipartite entanglement 
monotone that, as we will soon show, can be regarded as the continuous-variable tangle E T . 
For a generic pure state \ip) of a (1 + JV)-mode CV system, we define the square of the 
logarithmic negativity: 

This is a proper measure of bipartite entanglement, being a convex, increasing function of the 
logarithmic negativity Ejy, which is equivalent to the entropy of entanglement for arbitrary 
pure states. For any pure multimode Gaussian state \ip), with CM <r p , of 7V+1 modes assigned 
in a generic bipartition 1 x JV, explicit evaluation gives immediately that E t (i/j) = E T (cr p ) 
takes the form 

Er(o*) = In 2 fl/fi, - y/l/fi - lj , (10) 

where p\ = l/\/Det ct\ is the local purity of the reduced state of mode 1 with CM (T\. 
Def. (|9j is naturally extended to generic mixed states p of (N + l)-mode CV systems through 
the convex-roof formalism (see also Ref. 1221 where a similar measure, the convex-roof 
extended negativity, is studied). Namely, we can introduce the quantity 

E T (p) = inf V ft £ T (^), (11) 

{pi,i>i} ^ 

where the infimum is taken over all convex decompositions of p in terms of pure states {\ipi) }, 
and if the index i is continuous, the sum in Eq. il It is replaced by an integral, and the 
probabilities {p^} by a probability distribution ir(ip). Let us now recall that, for two qubits, 
the tangle can be defined as the convex roof of the squared negativity [22:1 (the latter being 
equal to the concurrence 1201 for pure two-qubit states). Here, Eq. il It states that the convex 
roof of the squared logarithmic negativity properly defines the continuous-variable tangle, 
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or, in short, the contangle E T (p), in which the logarithm takes into account for the infinite 
dimensionality of the underlying Hilbert space. 

From now on, we will restrict our attention to Gaussian states. For any multimode, 
mixed Gaussian states with CM er, we will then denote the contangle by E T (cr), in analogy 
with the notation used for the contangle E r (cr p ) of pure Gaussian states in Eq. dlQi . Any 
multimode mixed Gaussian state with CM er, admits at least one decomposition in terms of 
pure Gaussian states er p only. The infimum of the average contangle, taken over all pure 
Gaussian state decompositions, defines then the Gaussian contangle G T : 

G T (er) ee inf / n(d(T p )E T (<T p ) . (12) 

{TT(d<TP),<TP} J 

It follows from the convex roof construction that the Gaussian contangle G T (er) is an upper 
bound to the true contangle E T (er) (as the latter can be in principle minimized over a non- 
Gaussian decomposition): 

E T (a) < G r (cr) , (13) 

and it can be shown that G r (er) is a bipartite entanglement monotone under Gaussian local 
operations and classical communication (GLOCC) 1123 1 1241 . In fact, for Gaussian states, the 
Gaussian contangle, similarly to the Gaussian EoF 1231 . takes the simple form 

G T (er) = inf E t {<t p ) , (14) 

er p <er 

where the infimum runs over all pure Gaussian states with CM cr p < er. Let us remark that, 
if er denotes a mixed symmetric two-mode Gaussian state, then the Gaussian decomposition 
is the optimal one 1131 (it is currently an open question whether this is true for all Gaussian 
states 1251 ). and the optimal pure-state CM er p minimizing G T (er) is characterized by having 
n_ (er p ) = n_ (er) 11231 . The fact that the smallest symplectic eigenvalue is the same for both 
partially transposed CMs entails that E T (cr) = G r (er) = [max{0, — hin_(er)}] 2 . We thus 
consistently retrieve for the contangle, in this specific case, the expression previously found 
for the mixed symmetric reductions of fully symmetric three-mode pure states, Eq. l|8}. 



3. Monogamy inequalities and residual multipartite entanglement 

3.1. Monogamy inequality for all three-mode Gaussian states 

We are now in the position to prove the first main result of the present paper: all three-mode 
Gaussian states satisfy the monogamy inequality (|6j, using the Gaussian contangle G T (or 
even the true contangle E T for pure states) to quantify bipartite entanglement. 

We start by considering pure Gaussian states er p of three modes, each of the three reduced 
single-mode states being described respectively by the CMs er.;, crj, er fc . Due to the equality 
of the symplectic spectra across a bipartite cut, following from the Schmidt decomposition 
operated at the CM level |26|, any one of the two-mode reduced CMs cry, er^, (jjk, will 
have smallest symplectic eigenvalue of the associated partially transposed CM equal to 1, 
and will thus represent a mixed state of partial minimum uncertainty 1271 1281 . Given the three 
modes i, j, k, any such two-mode reduction, say that of modes i and j, is completely specified, 
up to local single-mode unitary operations, by the two local purities fii = (Det er,) -1 / 2 and 
/ij = (Det er,) -1 / 2 , and by the global two-mode purity equal to the local purity of mode 
ft: = (Dot er.y)" 1 / 2 = (Det er^.) -1 / 2 = fi k . Moreover, because the contangle between 
mode i and modes (jk) is a function of ^ alone, see Eq. ( II Ot . all the entanglement properties 
of three-mode pure Gaussian states are completely determined by the three local purities /x, , 
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Hj, and fXf., or by the associated local mixednesses Oj = 1/ Hi, cij = 1/fJ.j, and = l//ifc. 
The local mixednesses <Zj, a,-, and have then to vary constrained by the triangle inequality 

\a>j — Ofe| + 1 < a% < a,j +a k - 1 , (15) 

in order for cr p to be a physical state. This is a straightforward consequence of the uncertainty 
relation Ineq. Q applied to the reduced states of any two modes (see 1291 1241 for further 
details). Notice that Ineq. dl5> is a stronger requirement than the general Araki-Lieb inequality 
holding for the Von Neumann entropies of each single-mode reduced state 1129 1 . For ease of 
notation, let us rename the mode indices: {£, j,k} = {1, 2, 3}. Without loss of generality, we 
can assume a x > 1 (if a 1 = 1 the first mode is not correlated with the other two and all terms 
in Ineq. (|6} are trivially zero). Moreover, we can restrict to the case of er 12 and er 13 being both 

1 1 (23) 

entangled. In fact, if e.g. <Xi 3 denotes a separable state, then E T (a\<i) < E T \cr p ) because 
tracing out mode 3 is a LOCC (see Ref. II 11 ), and thus the sharing inequality is automatically 
satisfied. 

We will now prove Ineq. l|6) in general by using the Gaussian contangle, as it will imply 
the inequality for the true contangle; in fact Gr l(23) (cr' p ) = £^ l(23) (erf) but G^\a) > 
Er^{tr), I = 2, 3. Our strategy will be to show that, at fixed oi, i.e. at fixed entanglement 
between mode 1 and the remaining modes: 



= ln 2 ( ai - - 1) , (16) 

the maximum value of the sum of the 1|2 and 1|3 bipartite entanglements can never exceed 

Er^ 23 \ Namely, max Sl ^ Q < ln 2 (ai — \J a\ — 1), where 

Q = G r (<x 12 ) + G r (<T 13 ) , (17) 

and the variables a 2 and a 3 have been replaced by the variables s\ = (a 2 + a 3 )/2 and 
di = (oa — a 3)/2- The latter are constrained to vary in the region 

defined by the triangle inequality d!5l > and by the condition of the reduced two-mode 
bipartitions being entangled 1241 . We recall now that each an, I = 2, 3, is a two-mode state of 
partial minimum uncertainty. For this class of states the Gaussian measures of entanglement, 
including G T , can all be determined and have been computed explicitely 1241 . Skipping 
straightforward but tedious calculational details, and omitting from now on the subscript 1, 
we have that 

Q = In 2 [m(a, s, d) — y/ m 2 (a, s, d) — 1] 



+ In 2 [m(a, s, -d) - y / m 2 {a,s,d) - 1] , (19) 
where m = m_ if D < 0, m = m + otherwise, and 

m_= |fc_|/[( S -d) 2 -l] , 

m+ = {2 [2a 2 (l + 2s 2 + 2d 2 ) - (4s 2 - l)(4rf 2 - 1) - a 4 - VS]} 1/2 /[A{s - d)] , 



D = 2{s -d)- ^2[k 2 _+ 2k+ + \k-\(k 2 _+ 8k+ ) V2] /k + , 
k± = a 2 ±(s + d) 2 , 

5 = (a-2d-l)(o-2d+l)(a+2d-l)(a+2d+l)(a-2s-l)(o-2s+l)(a+2s-l)(a+2s+l) 
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Studying the derivative of with respect to s, it is analytically proven that, in the 
whole space of parameters {a, s, d} given by Eq. dl 81 . both m_ and m + are monotonically 
decreasing functions of s. The quantity Q is then maximized over s for the limiting value 
s = s mm = (a + l)/2. This value corresponds to three-mode pure states in which the 
reduced partition 2|3 is always separable, as it is intuitive because the bipartite entanglement 
is maximally concentrated in the 1|2 and 1 1 3 partitions. With the position s = ,s mm , D can be 
easily shown to be always negative, so that, for both reduced CMs <j\2 and <Ti3, the Gaussian 
contangle is defined in terms of m_ . The latter, in turn, acquires the simple form 

min n 1 + 3a + 2d 

m-{a,s ,d ) = — — — . 

6 + a — 2d 

Consequently, the quantity Q is immediately seen to be an even and convex function of d, 
which entails that it is globally maximized at the boundary \d\ = d max = (a — 1) /2. It turns 
out that 



Q[a, s = s min , d = ±d max ] = In 2 (a - V a 2 - 1) , (20) 

which implies that in this case the sharing inequality (|6j is saturated and the genuine tripartite 
entanglement is exactly zero. In fact this case yields states with = a\ and a.3 = 1 (if 
d = <i max ), or (13 = 01 and a-2 = 1 (if d = — d max ), i.e. tensor products of a two-mode 
squeezed state and a single-mode uncorrected vacuum. Being the above quantity Eq. (I20> the 
global maximum of Q, Ineq. holds true for any Q, that is for any pure three-mode Gaussian 
state, choosing either the Gaussian contangle G T or the true contangle E T as measures of 
bipartite entanglement. 

By convex roof construction, the above proof immediately extends to all mixed 
three-mode Gaussian states er, with the bipartite entanglement measured by G T . Let 
{Tr(dcrf n ), erf n } be the ensemble of pure Gaussian states minimizing the Gaussian convex 
roof in Eq. Jl 2b : then, we have 



> J 7r(d<)[G^«) + Gf«)] 



>Gf{<r) + Gf {a) , (21) 

where we have exploited the fact that the Gaussian contangle is convex by construction. This 
concludes the proof of the monogamy inequality for all three-mode Gaussian states. ■ 



3.2. Residual tripartite entanglement and monotonicity 

The sharing constraint © leads naturally to the definition of the residual contangle as a 
quantifier of genuine tripartite entanglement (arravogliamenf). This is in complete analogy 
with the case of qubit systems, except that, at variance with the three-qubit case (where the 
residual tangle of pure states is invariant under qubit permutations), for CV systems of three 
modes the residual contangle is partition-dependent according to the choice of the reference 
mode (but for the fully symmetric case). Then, the bona fide quantification of tripartite 
entanglement is provided by the minimum residual contangle: 



Et^ k = min 



Ei\W_ E i\i_ E i 



(22) 
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where (i,j, k) denotes all the possible permutations of the three mode indexes. This definition 
ensures that is invariant under mode permutations and is thus a genuine three-way 

property of any three-mode Gaussian state. One can verify that 

£*!(**) _ E l J k - (E^ - Ep k ) > 

if and only if a, > a,j, and therefore the absolute minimum in Eq. i22\ is attained by the 
decomposition realized with respect to the reference mode i of smallest local mixedness dj, 
i.e. of largest local purity /ij. 

A crucial requirement for Er to be a proper measure of tripartite entanglement is 
that it be nonincreasing under LOCC. The monotonicity of the residual tangle was proven 
for three-qubit pure states in Ref. 1301 . In the CV setting, we will now prove that for pure 
three-mode Gaussian states the residual Gaussian contangle G$^ k , defined in analogy with 
Eq. \22\ . is an entanglement monotone under tripartite GLOCC, and specifically that it is 
nonincreasing even for probabilistic operations, which is a stronger property than being only 
monotone on average 11231 . We thus want to prove that 

G^ k (G p ((r p )) < Gf lk (cr p ) 

, where G p is a pure GLOCC mapping pure Gaussian states cr p into pure Gaussian states 
13 II 1321 . Every GLOCC protocol can be realized through a local operation on one party only. 
Assume that the minimum in Eq. (1221 is realized for the reference mode i\ the output of a 
pure GLOCC G p acting on mode i yields a pure-state CM with a[ < ai, while cij and at 
remain unchanged 113 II . Then, the monotonicity of the residual Gaussian contangle Gr^ k 
under GLOCC is equivalent to proving that G % }^ k = Gt^^ — Gr J — G % } k is a monotonically 
increasing function of for pure Gaussian states. One can indeed show that the first derivative 
of Gr^ k with respect to a^, under the further constraint aj < a^, is globally minimized for 
a,i = cij = <2fc = a; oc , i.e. for a fully symmetric state. It is easy to verify that this minimum 
is always positive for any ai oc > 1, because in fully symmetric states the residual contangle 
is an increasing function of ai oc . Therefore the monotonicity of G 1 }^ under GLOCC for all 
pure three-mode Gaussian states is finally proven. ■ 



3.3. Monogamy inequality for N-mode symmetric Gaussian states 

We next want to investigate whether the monogamy inequality can be generalized to 
Gaussian states with an arbitrary number N + 1 of modes, namely whether 



N 



E 



*ICj'i 



- E%bi > ° 



(23) 



i=i 



Establishing this result in general is a highly nontrivial task, but we can readily prove it for all 
symmetric multimode Gaussian states. As usual, due to the convexity of G T , it will suffice to 
prove it for pure states, for which the Gaussian contangle coincides with the true contangle in 
every bipartition. For any N and for a; oc > 1 (for a; oc = 1 we have a product state), 

Ei\Ui,-M = ln 2 (a;oc _ ^a 2 loc -l) (24) 

is independent of N, while the total two-mode contangle 



NE^' = — In 2 
4 



3 (JV + 1) - 1 



(«L-i)(«L(^ + i) 2 -(^-i) 2 ) In 



(25) 
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Figure 1. Behaviour of the Gaussian contangle G T between mode 1 and all the other 
remaining three modes (horizontal axis), plotted versus the total bipartite Gaussian contangle 
2j=2 (vertical axis). The two quantities have been evaluated numerically in 60000, 
randomly generated, four-mode Gaussian states. States saturating the monogamy inequality 
Eq. 1231 fall on the solid red line. No physical state has been found lying in the upper-left 
half of the plane, above the red line. This numerical result supports the conjecture that the 
Gaussian contangle is monogamous, and satisfies Eq. 1231 in all multimode Gaussian states. It 
is important to notice that, for increasing entanglement between one mode and the others, the 
states tend to distribute farther away from the 45-degree line boundary, hinting at the presence 
of genuine multipartite entanglement shared between all the modes. 

is a monotonically decreasing function of the integer N at fixed a\ oc . Because the sharing 
inequality trivially holds for N = 1, it is inductively proven for any N. ■ 
This result, together with extensive numerical evidence obtained for randomly generated 
non-symmetric 4-mode Gaussian states (see Fig.Q, strongly supports the conjecture that the 
monogamy inequality be true for all multimode Gaussian state, using the (Gaussian) contangle 
as a measure of bipartite entanglement . It is currently under way to fully prove the conjecture 
analytically. 

4. Structure of entanglement sharing: the CV GHZ/W states 

We are now in the position to analyze the sharing structure of CV entanglement by taking 
the residual contangle as a measure of tripartite entanglement, in analogy with the study 
of three-qubit entanglement 1301 . Namely, we pose the problem of identifying the three- 
mode analogues of the two inequivalent classes of fully inseparable three-qubit states, 
the GHZ state |33) |V>ghz) = (l/v^) [|000) + |111)], and the W state M \ipw) = 
(1/V3) [|001) + |010) + 1 100)] . These states are pure and fully symmetric. On the one hand, 
the GHZ state possesses maximal tripartite entanglement, quantified by the residual tangle 
I3l l30l . without any two-qubit entanglement. On the other hand, the W state contains the 
maximal two-party entanglement between any pair of qubits 1 30| and its tripartite residual 
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tangle is consequently zero. 

Surprisingly enough, in symmetric three-mode Gaussian states, if one aims at 
maximizing (at given single-mode squeezing r; oc or single-mode mixedness a; oc ) either the 
two-mode contangle in any reduced state {i.e. aiming at finding the CV analogue of the 
W state), or the genuine tripartite contangle (i.e. aiming at defining the CV analogue of the 
GHZ state), one finds the same, unique family of pure symmetric three-mode squeezed states. 
These states, previously named "GHZ-type" states 1191 . can be defined for generic iV-mode 
systems, and constitute an ideal test-ground for the study of multimode CV entanglement 
1 1 81 1 1 01 . The peculiar nature of entanglement sharing in this class of CV GHZ/W states is 
further confirmed by the following observation. If one requires maximization of the 1x2 
bipartite contangle E T under the constraint of separability of all two-mode reductions, 
one finds a class of symmetric mixed states whose tripartite residual contangle is strictly 
smaller than the one of the GHZ/W states, at fixed local squeezing. Therefore, in symmetric 
three-mode Gaussian states, when there is no two-mode entanglement, the three-party one is 
not enhanced, but frustrated. Instead, if a mode is maximally entangled with another, it can 
also achieve maximal quantum correlations in a three-party relation. 

These results, unveiling a major difference between discrete-variable and CV systems, 
establish the promiscuous nature of CV entanglement sharing in symmetric Gaussian states. 
Being associated to degrees of freedom with continuous spectra, states of CV systems need 
not saturate the CKW inequality to achieve maximum couplewise correlations. In fact, 
without violating the monogamy inequality (|6jl, pure symmetric three-mode Gaussian states 
are maximally three-way entangled and, at the same time, maximally robust against the loss 
of one of the modes due, for instance, to decoherence. 

Furthermore, the residual contangle Eq. J22b in GHZ/W states acquires a clear operative 
interpretation in terms of the optimal fidelity in a three-party CV teleportation network 1341 . 
This finding readily provides an experimental test, in terms of success of teleportation- 
network experiments I35II36I . to observe the promiscuous distribution of CV entanglement in 
symmetric, three-mode Gaussian states. 

5. Concluding remarks 

It is a central trait of quantum information theory that there exist limitations to the free sharing 
of quantum correlations among multiple parties |4|. This aspect can be quantified in terms of 
monogamy constraints, as first introduced by Coffman, Kundu and Wootters for states of three 
qubits. In this paper we have generalized these monogamy constraints to infinite-dimensional 
systems. This extension required the definition of a proper entanglement monotone, able to 
capture the trade-off between the couplewise entanglement and the genuine tripartite and, in 
general, multipartite entanglement in multimode Gaussian states. We proved analytically that 
the continuous-variable entanglement is monogamous in all three-mode and in all symmetric 
multimode Gaussian states, and have numerically convincing evidence that this holds true in 
all TV-mode Gaussian states as well. 

Very remarkably, in the case of pure states of three modes, the residual entanglement 
emerging from the monogamy inequality turns out to be a tripartite entanglement monotone 
under Gaussian LOCC, representing the first bona fide measure of genuine multipartite 
entanglement for CV systems. This measure has been applied to investigate the concrete 
structure of the distributed entanglement in three-mode Gaussian states, leading to the 
discovery that there exists a special class of states (pure, symmetric, three-mode squeezed 
states) which simultaneously maximize the genuine tripartite entanglement and the bipartite 
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entanglement in the reduced states of any pair of modes. This property, which has no 
counterpart in finite-dimensional systems, has been labeled by us promiscuous sharing of 
CV entanglement. 

The collection of results presented here is of basic importance for the understanding 
of quantum correlations among multiple parties in systems with infinitely many degrees of 
freedom. Several hints emerging from our results are worth being investigated. Among 
them, the analysis of the effect of decoherence on states with promiscuous entanglement 
sharing, to study the actual robustness of this correlation structure in the presence of noise 
and mixedness. From a more fundamental point of view, the extension of the present 
techniques to define a more general measure of multipartite entanglement for all multimode, 
Gaussian and non-Gaussian states is the most evident, although very challenging, task to 
aim at. Another important follow-up should concern the experimental implications of our 
findings, ranging from the verification of the promiscuity in terms of quantum communication 
experiments 1341 . to the preparation of special classes of multiparty entangled Gaussian states 
and the implementation of these resources for novel protocols of quantum information with 
continuous variables. All these issues are being currently and actively studied 1291 . 

In a wider perspective, we wish to mention that monogamy inequalities should play 
a fundamental role in the understanding of entanglement and other structural properties 
of a wide class of complex quantum systems both in discrete and continuous variables 
1161 1371 l38l 1391 Hoi . At this stage, in the continuous-variable setting, the most interesting 
and urgent open problems are perhaps the extension of the CKW-type sharing inequalities to 
non-Gaussian states and to systems in which the reference party is constituted by more than 
one mode. The latter problem is the CV analogue of extending the monogamy inequalities 
to systems in which the reference party is constituted by more than one qubit 1161 . Finally, 
another very intriguing open issue is the generalization to systems of arbitrary finite dimension 
D interpolating between the case D = 2 and the infinite-dimensional instance (4j. In all 
these problems the challenge is to define easily computable measures of entanglement that 
allow concrete quantifications and bounds on the sharing structure and distribution of quantum 
correlations. 
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